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ABSTRACT 



As part of the long-term radial velocity monitoring of known transiting planets - 
designed to measure orbital eccentricities, spin-orbit alignments and further planetary 
companions - we have acquired radial velocity data for the two transiting systems 
WASP-12 and WASP-14, each harbouring gas giants on close orbits (orbital period of 
1.09 and 2.24 days respectively). In both cases, the initial orbital solution suggested 
a significant orbital eccentricity, 0.049 ± 0.015 for WASP-12 and 0.091 ± 0.003 for 
WASP-14. Since then, measurements of the occultation of WASP-12 in the infrared 
have indicated that one projection of the eccentricity (e cosw) was close to zero, casting 
doubt on the eccentricity from the initial radial velocity orbit. Our measurements 
confirm that the initial eccentricity detection could be spurious, and show that the 
radial velocity data is compatible with a circular orbit. A MCMC analysis taking into 
account the presence of correlated systematic noise in both the radial velocity and 
photometric data gives e = 0.017^Q gi[^. By contrast, we confirm the orbital eccentricity 
of WASP-14, and refine its value to e = 0.088 ± 0.003. WASP-14 is thus the closest 
presently known planet with a confirmed eccentric orbit. 
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1 INTRODUCTION 

Transiting planets are an important source of information on 
the formation, structure and evolution of extra-solar planets. 
We are monitoring known transiting planetary systems in ra- 
dial velocity with the SOPHIE spectrograph in the Northern 
hemisphere and HARPS spectrograph in the Southern hemi- 
sphere to refine our knowledge of the dynamics of these sys- 
tems, notably the orbital eccentricity, spin-orbit angle and 



* Based on observations made at the 1.93-m telescopes at Ob- 
servatoire de Haute-Provence (CNRS), Prance with the SOPHIE 
spectrograph. 



presence of additional companions (e.g. [Loeillet et al.|2008| 
Hebrard et aLpOOS, ESO Prog. 0812.C-0312) 



In this paper, we analyse our new SOPHIE radial- 
velocity data for two transiting planetary systems, WASP- 
12 and WASP-14. Both are characterized by close-in but 
apparently eccentric orbits, and therefore represent poten- 
tially important systems to constrain the migration, tidal 
and thermal evolution of gas giant planets. We combine our 
radial-velocity data with previously published data and a 
realistic treatment of correlated noise to calculate updated 
constraints on the orbital eccentricities. 

The companion of the 11.7th-magnitude star WASP- 
12 is a particularly interesting example ( |Hebb et al.||20d9} 
hereafter H09). It orbits extremely close to its host star. 
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even by the standards of the so-called "hot Jupiters" , with 
a period of 1.09 days, corresponding to an orbital distance 
only 3 times the radius of its host star. Moreover, WASP-12b 
has an inflated radius, R~ 1.8 Rj, one of the most extreme 
examples of anomalous radii for hot Jupiters. As a result, 
the planet fills about half of its Roche lobe ( Li et al.|2010 |. 

With such a short orbital distance and large size, a gas 
giant planet is expected to undergo complete orbital syn- 
chronisation and circularisation on a short timescale, much 
shorter than the age of a typical field main-sequence cool 
star. Indeed, most planets orbiting closer than 0.05 AU are 
observed to have circular orbits. However, H09 determined 
a value of e = 0.049 ± 0.015 for the orbital eccentricity 
of WASP-12, a significant departure from circularity. This 
would make the planet by far the subject of the strongest 
tidal dissipation in any known planetary system. The mea- 
sured eccentricity is based on fitting a Keplerian orbital mo- 
tion on the radial velocity measurements collected by H09 
with the SOPHIE spectrometer ( [Perruchot et al.||2008| ) to- 
gether with transit photometry. 



Li et al. (20101 studied the case of WASP-12 with that 



value of eccentricity, and found a large implied mass loss 
and dissipation of tidal energy in the planet. 

In a transiting system, the time lag between the tran- 
sit and the occultation has a strong dependence on the 
projected orbital eccentricity (ecoso;). Therefore, if the oc- 
cultation can be detected with sufficient significance, this 
provides a stringent test of the eccentricity. |Lopez-Morales| 
et al. (2009 hereafter L09) have measured the occultation 
of WASP-12b from the ground with SPICam on the ARC 
telescope at Apache Point Observatory in the z' band. Their 
best-fit result indicated a occultation with a significant time 
lag compared to the epoch expected for a circular orbit, with 
a similar level of significance to H09. Nevertheless, the pres- 
ence of residual correlated noise is apparent in the L09 data 
(see Fig. [2|, as expected for ground-based photometry at 
such a high accuracy - the depth of the occultation is only 
about 0.08±0.02 %. 

As a result, the issue remained inconclusive until a 
space-based measurement of the occultation with the Spitzer 
Space Telescope ( jCampo et al.|[20To| hereafter CIO) unam- 
biguously showed that the timing of the occultation was 
precisely that expected for a circular orbit. This result sug- 
gested that the L09 time lag was probably due to instrumen- 
tal systematics, and that the orbit of WASP-12 was probably 
circular, since a fine-tuned alignment would be required to 
reconcile the Spitzer result with the H09 value of the eccen- 
tricity. 

It is interesting to note that there is an inherent bias 
in eccentricity measurements from radial velocities, because 
a Keplerian orbit cannot get more circular than e = 0. 
Any noise applied to a circular orbit will result in an ec- 
centric best-fit orbit. Underestimating the noise will lead 
to spurious detections of small eccentricities. This was al- 
ready recognized in the context of stellar binaries by |Lucy| 
& Sweeney (19711. These authors showed that spurious ec- 



centricity detections tended to dominate for e < 0.1 for a 
typical precision at that time and stellar binary amplitudes. 
Four decades later, both companion masses and RV accu- 
racies having changed by about three orders of magnitudes, 
and the same issue resurfaces for exoplanets. 

WASP-14 is, after WASP-12, the known transiting 



planet having a reported non-circular orbit (Joshi et al 
[2009, e = 0.091 ± 0.003) with the second-shortest period 
(P=2.2 days). This makes it another test-case for tidal evolu- 
tion of close-in gas giants. If its orbital eccentricity is indeed 
near 0.1, then this non-zero but relatively low value - in the 
context of the distribution of giant exoplanet eccentricities 
- makes it likely that this planet has undergone some degree 
of orbital evolution, and is still subject to strong tidal forces 
at present. Therefore its presence may be useful to constrain 
the tidal synchronisation timescale. It is also an important 
object when studying the issue of the anomalous radius of 
hot Jupiters because of its infiated size, with Rp — 1.28 Rj. 
WASP-14 occupies a distinctive position in the relevant pa- 
rameter space: irradiation, orbital distance, eccentricity and 
size. 



2 OBSERVATIONS 

We obtained 29 radial- velocity measurements for WASP-12 
(16 during a single night, and 13 at various values of orbital 
phase) and 11 for WASP-14, using the SOPHIE spectro- 
graph installed on the 1.93-m telescope at OHP (France). 
The observations were gathered between 17 January 2009 
and 27 March 2010. The 16 in-transit measurements for 
WASP-12 were obtained with the objective of constraining 
the spin-orbit angle via the Rossiter-McLaughlin effect. 

SOPHIE is a spectrograph optimized for precise radial- 
velocity measurements and has participated in the detec- 
tion of numerous transiting exoplanets in the northern 
hemisphere, notably from the WASP and CoRoT transit 
searches. It reaches a stability of a few ms~^ for bright tar- 
gets. WASP-12 and WASP-14, however, are near the faint 
end of the capacity of the 1.93-m telescope, and were mea- 



sured in the "High Efficiency" mode of SOPHIE (See Perru- 
chot et al.[ 2008 Bo uchy et al.|2009 1. This mode has a higher 
throughput than the standard mode, the "High Resolution" 
mode, thus allowing fainter targets to be measured, but is 
less optimized for radial velocity. When considering the en- 
semble of data for known transiting planets obtained with 
SOPHIE, we have found evidence for large excursions of the 
velocity zero-point with time (to the level of several dozen 
ms~^ in some cases). As part of the constant improvement 
of the SOPHIE reduction pipeline, this effect is monitored 
and corrected for as far as possible, but the presence of rela- 
tively large instrumental systematics in the High Efficiency 
data is a possibility, especially with older data collected be- 
fore we became aware of the issue. This must be remembered 
when performing an orbital analysis based on data from the 
High Efficiency mode. 



3 ANALYSIS AND RESULTS 

The orbital parameters of the two transiting planets were 
calculated from the radial velocity data, (together with 
published photometry data for the transit and occulta- 
tion in the case of WASP-12), with a Marko Chain Monte 
Carlo (MCMC) method. The main advantage of the MCMC 
method is that it allows a seamless combination of radial- 
velocity data with light-curve data both for the transit and 
occultation, as well as information on the parent star. The 
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use of MCMC in this context is described by, among oth- 
ers, Holman et al. (20061. Our implementation is described 



H09 transit to Ua ~ 0.1274 and u^, — 0.3735 according to 



Claret (20041. 



Pont et al. (20091. We use a Bayesian treatment of the 



a priori constraints from stellar evolution models with the 



method of Pont & Eyer (20041. We model the radial veloc- 



ity using a Keplerian orbit, the Rossiter-McLaughlin effect 



on the radial velocity using Gimenez (20061 and the transit 



and eclipse light curves using Mandel & Agol ( 2002 1 



MCMC methods are powerful, but they also tend to 
obscure the relation between the errors on the measured 
data and the calculated uncertainties on the final system 
parameters. Neglecting non-random sources of noise in the 
data (such as instrumental systematics and stellar variabil- 
ity) can lead to an underestimation of the uncertainties in 



the final system parameters by a large factor (Pont et al 



2006 1, and, in the case of orbital eccentricities, to a system- 
atic bias ( [Lucy fc Sweeney||1971[ ). We account for the pres- 
ence of correlated noise in both the photometric and ra- 
dial velocity data by modifying the merit function used in 
the MCMC to include the possible presence of non-random 
noise. 



3.1 The orbital eccentricity of WASP-12 

Figure [l] displays our radial-velocity data for WASP-12, to- 
gether with a circular orbital solution (solid line, residuals 
shown in the lower panel) and the eccentric-orbit fit of E[09 
(dotted). Our radial velocity data is shown together with 
the H09 data in the right panel. As the Figure shows clearly, 
the radial-velocity signal cannot be adequately modelled by 
a periodic orbital signal affected by random noise. The pres- 
ence of a correlated non-periodic component is especially 
obvious during the in-transit sequence. This could be due 
to instrumental noise, stellar variability or an unaccounted 
planetary companion in the system, all of which would be- 
have in the same way as far as the orbital fit is concerned. 
Given our experience with SOPHIE in High Efficiency mode, 
we consider the first cause as likely. 

We use the transit lightcurve from H09 and the eclipse 
lightcurve from L09 with our own radial velocity measure- 
ments to fit a Keplerian orbit, and we account for the effect 
of red noise by modifying the uncertainties on the data us- 
ing cr' = [a^ -\- Na^)^^^, where a is the random error, a' 
the modified error used to compute the merit function for 
the MCMC, CTr the red- noise factor, and the number of 
data points over a typical correlation timespan. The reader 
is referred to Pont et al.| (|2006) and Winn (2010 in prep) 
for more details on this approach. For the H09 photometric 
transit, we estimate (Tr = 0.0005 and A'^ = 20. Similarly, for 
the L09 data, we set ar = 0.0002 and A = 9. For our own 
radial velocity data, we estimate the red noise parameter to 



be 8 ms~ , and A' = 16 during the transit and A = 1 out- 
side the transit (these points are taken in different nights so 
that the errors are not expected to be correlated). 

We vary the period P, mid-transit time T^j., system ve- 
locity vo, eccentricity components ecostj and esinoj, semi- 
amplitude K, impact parameter 6, scaled radius Rp/Ra, the 
mass Ms, radius Rs and T^f^ '^^ the star. The scaled semi- 
major axis a/Ra is calculated from the assumption of a Ke- 
plerian orbit using the period P. 

We set the quadratic limb darkening parameters for the 



Our results for the orbital parameters of the WASP-12 
system are shown in Table [l] Table [l] shows the results of 
our new SOPHIE RV data used in combination with the H09 
transit lightcurve, and the L09 and CIO eclipse lightcurves, 
compared to the H09 parameters. 

Although our data covers a complete spectroscopic 
transit, and thus potentially constrains the projected spin- 
orbit angle of the WASP-12 system through the Rossiter- 
McLaughlin effect, the constraint is weak once the possible 
presence of non-random noise is taken into account. The dis- 
tribution of the spin-orbit angle from our MCMC spans a 
wide interval extending from a prograde orbit to a projected 
spin-orbit angle larger than 90 degrees. The data marginally 
favours a prograde rather than retrograde orbit. 

The best solution with our radial velocity data and the 
constrains from the transit and eclipse light curves gives 
e — 0.017^0010- This is a marginal detection, so we set 
e = to work out the rest of the system parameters. The H09 
radial-velocity data forces the solution towards a higher ec- 
centricity, significant at the ~3cr level, but Figure [l] strongly 
suggests that this is probably an artefact due to an excursion 
of the zero-point between different nights. 

It is interesting to compare the situation in radial ve- 
locity with the similar sequence of events regarding the oc- 
cultation photometric data. Figure [5] (left panel) shows the 
eclipse data for both L09 (top two panels) and CIO (bottom 
two panels). The solid line represents the L09 solution, and 
the residuals are plotted for both datasets. The L09 solu- 
tion is somewhat plausible for the L09 data but definitely 
not for the CIO data. The dotted line shows the a circular 
orbit e = 0. Figure [5] (right panel) shows the same data, but 
the solid line is now a circular orbit e = 0, and the residuals 
for both datasets are again plotted. It is clear that the cir- 
cular solution fits the CIO dataset and remains reasonable 
for the L09 dataset. This could be explained if the effects of 
instrumental systematics had been underestimated in L09. 

3.2 The orbital eccentricity of WASP-14 

Figure[3|s hows our SOPHIE data and the FIES and SOPHIE 
data from jJoshi et al.| ( |2009[ ) along with with the best-fit or- 
bit and a circular orbit. We adopt the prior distribution on 



the period P from photometric data by Johnson et al. ( 2009 \ 
and that on the mid-transit time Ttr from |Joshi et al.| ( |2009[ ). 
We include their value of Rp/Rs, a/ Rs and v sini to compute 
the Rositter-McLaughlin effect, and estimate the quadratic 
limb darkening coefficients Ua = 0.405 and Ub = 0.337 from 



Claret (20041. We work out the orbital parameters using 



both our radial velocity data and that published in the dis- 
covery paper (since we do not have any radial velocity data 
during the transit, the RM effect from that paper are not 
significantly modified). We include the effect of red noise in 
a similar manner explained above, adding 8 ms~^ of cor- 
related uncertainties to each single datapoint and \/N x 8 
ms"'^ to each point in each group of N points taken on the 
same night. 

Our results are shown in Table (2] Our best-fit value 
for the orbital eccentricity is e = 0.088 ± 0.003, in good 
agreement with the value found by |Joshi et al.| ( |2009[ ) . 

We examined the possibility of a scenario similar to 
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Figure 1. Left— Plot showing our new SOPHIE radial velocity data for WASP-12, overplotted with a circular orbit (solid line) and an 
orbit with the eccentricity in H09. The residuals relative to the circular orbit are shown in the bottom panel. 

Right— Plot showing all the radial velocity data for WASP-12, with our new SOPHIE data plotted with square symbols and the H09 
data plotted with dots. Residuals for an orbit with the H09 value of eccentricity, e = 0.049, are shown in the bottom panel. 



Parameter H09 New SOPHIE RV and photometry 



Centre-of-mass velocity Vq [km s ■'-] 


19.085±0.002 


19.061±0.007 


Orbital period P [days] 


1.091423±0.000003 


1.09139±0.00001 


Orbital eccentricity e 


0.049±0.015 


OQ^y-f 0.015* 

'J-'-'^' -0.010 


Argument of periastron lj [°] 


-74+" 


(unconstrained) 


Velocity semi-amplitude K [ms~^] 


226±4 


233±6 


Epoch of transit Ttr [BJD] 


2454508. 9761±0.0003 


2454508. 9761±0.0006 



Table 1. System parameters for WASP-12. Left: H09. Right: Our SOPHIE radial velocity data, H09 transit photometry data and L09 
and CIO eclipse photometry data. 

*We adopt a circular orbit for working out the rest of the system parameters. 



that of WASP-12, with correlated noise causing a spurious 
eccentricity detection. The lower panel of Figure [3] shows 
the residuals around the best-fit circular orbit (given by a 
MCMC run with e = fixed). We find that, in contrast to 
the case of WASP-12, the differences between observation 
and model assuming a circular solution are periodic and reg- 
ular, which would not be the case for correlated noise. We 
also looked for possible evidence for a second planetary com- 
panion in the system by examining the residuals as function 
of time, finding no unambiguous trend. 



4 DISCUSSION 

Our results confirm the strong indications of CIO that all 
the available data for WASP-12 is compatible with a circular 
orbit, and that the eccentricity of the best-fit orbit to the 
radial velocity of H09 and subsequently the occultation data 
of L09 may be due to correlated noise. Not accounting for 
this noise in the statistical analysis could lead to an apparent 
~ 3cr significance for the rejection of the null hypothesis 
(e = 0), but the new data strongly suggest that the orbit of 
WASP-12b is indeed circular. 



CIO suggested that, for an eccentric orbit, the difference 
between the eclipse phase in CIO and L09 could be due to 
apsidal precession. This would require that the argument 
of periastron had changed from to = —74:'^\q in Feb 2008 
(H09) to w -90° in Oct 2008 (CIO). In terms of the 
projected component of the eccentricity, ecosuj, one obtains 
ecoscj = 0.014 ± 0.004 in Feb 2008 using the H09 orbital 
parameters, CIO found e coscj = 0.0019±0.0007 in Oct 2008, 
L09 found e cos a; = 0.0156±0.0035 in Feb-Oct 2009, and our 
resuh is e cos a; = 0.0037± 0.0035 for January 2009 to March 
2010. When considered in the order that observations were 
made, these values do not support the hypothesis of apsidal 
precession. 

A circular orbit for WASP-12 removes the needs for 
models to explain the survival of such an eccentricity at 
this very short period, in face of what would have been ex- 
tremely strong tidal effects. In particular, the scenario of O] 



et al. (20101, using the eccentricity from H09 to infer val- 



ues of mass loss and tidal dissipation for WASP-12, lose its 
principal empirical support. 

The eccentricity of WASP-14b, by contrast, is solidly 
confirmed by our measurements. This illustrates the capac- 
ity of SOPHIE to measure accurate values of orbital eccen- 
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Figure 2. Eclipse flux for WASP-12. Left— Top two panels; L09 flux and residuals. Bottom two panels: CIO flux and residuals. Solid 
lines represent the eclipse prediction for an eccentric orbit (L09 et al. e = 0.057), and dotted lines represent the eclipse prediction for a 
circular orbit (e = 0). Residuals are plotted for the L09 solution of e = 0.057, which is clearly inappropriate for the CIO eclipse (bottom). 
Right— Top two panels: LOO flux and residuals. Bottom two panels: CIO flux and residuals. Solid lines represent the eclipse prediction 
for a circular orbit (e = 0) and dotted lines represent the eclipse prediction for an eccentric orbit (LOO e = 0.057). Residuals are plotted 
for the circular orbit. 



Parameter 


Joshi et al 


This paper 


Centre-of-mass velocity Vb[km s^^] 




-4.985±0.003 


Orbital eccentricity e 


0.091±0.003 


0.088±0.003 


Argument of periastron u) [°] 


-106.6±0.7 


-107±1 


Semi-amplitude K [ms^^] 


993±3 


991±3 



Table 2. System parameters for WASP-14. (Left: Joshi et al. Right: Our results.) 



tricity for transiting planets, given a sufficient number of 
measurements well distributed in phase and spread over dif- 
ferent nights (the measurements for WASP-12 having most 
weight towards an eccentric solution were gathered during 
only two different nights). 



Eccentricity distribution and tidal circularisation 
for hot Jupiters 

Figure [4] shows a plot of the eccentricity of known transit- 
ing planets with a < 0.1 AU against their semi-major axis a 
(log scale). There is a well-know trend for the inner plan- 
ets to have circular orbits, due to tidal orbital decay. If the 
WASP-12 system follows a circular orbit, the WASP-14 sys- 
tem establishes a new lower limit for the point at which 
orbital eccentricity can survive tidal evolution for a suffi- 
ciently long time to be observed in a sample not selected 
by age. The fact that the eccentricity of WASP-14 is clearly 
non-zero, but also near the lower end of the eccentricity dis- 
tribution for planets out of reach of tidal evolution (broadly 
distributed between and 1) probably indicates that this ob- 
ject has already undergone some tidal evolution. As a result, 
WASP-14 is an interesting object to constrain the tidal evo- 
lution, intensity of tidal torque and relation between tides 
and size, for gas giant planets. 

WASP-10 is a third known transiting planets closer than 



0.04 AU with a possible non-circular orbit. However, the ec- 
centricity found by Christian et al. (2009) is significant to a 
level similar to H09 for WASP-12, and based in a comparable 
number and quality of RV measurements. The evidence for a 
circular orbit actually rests on two SOPHIE measurements, 
with residuals compared to a circular orbit that are compa- 
rable to the amplitude of the instrumental effects found in 
the High Efficiency mode. It is therefore possible that this 
detection too is due to underestimated instrumental system- 
atics acting on a circular orbit, in the manner described in 



Lucy & Sweeney ( 1971 1 



The fact that WASP-12b has a large size {Rp = 1.79 Rj, 
for a = 0.0229 and e ~ 0) and WASP-14b a smaller but still 
markedly inflated size (f?p = 1.28 Rj, for a = 0.037 AU and 
e ~ 0.09), provides a challenging test for theories attempting 
to account for the anomalous radii of some hot jupiters. 
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Figure 3. The new SOPHIE r adial velocity data for WASP-14 are 
shown with squares, while the'Joshi et al. I {20091) data from FIES 
and SOPHIE are shown with points and crosses respectively. The 
solid line is the best solution with e = 0.088, and the dotted line 
is a circular orbit. The middle panel shows residuals plotted for 
the best-fit orbit e=0.088, while the bottom panel shows residuals 
for a circular orbit ( (using the parameters of |Joshi et al.|2009l 
for the spectroscopic transit)). 
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Table 3. Radial velocity measurements for WASP-12 (errors in- 
clude random component only) 
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Figure 4. Plot of eccentricity e against semi-major axis a (log 
scale) for close-in transiting planets. Orbits with measured eccen- 
tricities are shown with a plus sign, and orbits with assumed zero 
eccentricities with a cross. 
* [Christian et al.| pOOOl l. 
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